Delay in atomic photoionization 
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We analyze the time delay between emission of photoelectrons from the outer valence ns and np 
sub-shells in noble gas atoms following absorption of an attosecond XUV pulse. By solving the time 
dependent Schrodinger equation and carefully examining the time evolution of the photoelectron 
wave packet, we establish the apparent "time zero" when the photoelectron leaves the atom. Various 
processes such as elastic scattering of the photoelectron on the parent ion and many-electron corre- 
lation affect the quantum phase of the dipole transition matrix element, the energy dependence of 
which defines the emission timing. This qualitatively explains the time delay between photoemission 
from the 2s and 2p sub-shells of Ne as determined experimentally by attosecond streaking [Science 
328, 1658 (2010)]. However, with our extensive numerical modeling, we were only able to account 
for less than a half of the measured time delay of 21 ± 5 as. We argue that the XUV pulse alone 
cannot produce such a large time delay and it is the streaking IR field that is most likely responsible 
for this effect. 

PACS numbers: 32.30.Rj, 32.70.-n, 32.80.Fb, 31.15.ve 



The attosecond streaking has made experimentally ac- 
cessible the characteristic timescale of electron motion in 
atoms [1, 2]. Among other spectacular applications of 
this technique, it has become possible to determine the 
time delay between subjecting an atom to a short laser 
pulse and subsequent emission of the photoelectron. In 
a recent work by Eckle et al [3], the helium atom was 
subjected to a near-infrared laser pulse with intensity of 
several units of 10 14 W/cm 2 . Such a strong field ion- 
ization regime could be characterized by a fairly small 
Keldysh parameter 7 ~ 1. The time delay in such a 
photoemission process can be conveniently analized in 
terms of nonadiabatic tunneling [4]. In a subsequent ex- 
periment by Schultze et al [5], the time delay was mea- 
sured in neon in the XUV photon energy range by high- 
order harmonic conversion of the driving near-infrared 
laser pulse. In this regime, which is characterized by a 
moderate intensity, short wavelength and 7 3> 1, it is 
believed that the formation of the outgoing wave packet 
follows instantaneously temporal variation of the incident 
electromagnetic field. Nevertheless, a sizable time delay 
of 21 ± 5 as was reported between photoemission from 
the 2s and 2p valence sub-shells of Ne. Schultze et al 
[5] argued that a comprehensive temporal characteriza- 
tion of photoemission on the attosecond time scale could 
provide a new insight into intra-atomic electron correla- 
tions. Indeed, the best theoretical treatment within an 
independent electron model could only account for 4.0 as 
time delay. When the theoretical model was corrected for 
electron correlations before and after photoionization, a 
relative delay of 6.4 as was obtained. This unresolved dif- 
ference between the measured and calculated time delays 
puts many-electron models of atomic photoionization un- 
der significant strain. If substantiated, this difference 
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could potentially point to new physical mechanisms un- 
derpinning electromagnetic interaction in atoms on the 
attosecond time scale. 

In this Letter, we perform extensive study of the time 
delay between the ns and np 2 outer valence sub-shell 
photoionization in noble gas atoms. We employ both the 
explicit time-dependent and stationary treatment of the 
photoionization process. To this end, we solve the time- 
dependent Schroding equation (TDSE) in the single ac- 
tive electron approximation. By carefully examining the 
time evolution of the photoelectron wave packet, we es- 
tablish the apparent "time zero" when the photoelectron 
leaves the atom. To account for electron correlation, we 
solve a set of coupled integral equations in the random 
phase approximation with exchange (RPAE) [6] . A sim- 
ilar approximation is employed to account for a possible 
target polarization effect [7]. 

Within an independent electron approximation, the 
time delay is caused by the energy dependence of the 
elastic scattering phase shifts of the photoelectron mov- 
ing in the Hartree-Fock (HF) potential of the Ne + ion. 
The many-electron correlation, which is due to inter-shell 
2s-2p coupling, depends on the energy of the photon and 
thus add additional component to the quantum phase of 
the dipole matrix element. Both these effects account for 
the time delay not exceeding ~ 10 as. This recovers only 
about one half of the experimental value of 21 ± 5 as. 
We carefully examine other correlation and polarization 
corrections, but find them unable to produce any sizable 
contribution to the measured time delay. We also analyze 
these effects in other noble gas atoms. 

The time-dependent calculation of photoionization in 
Ne was performed by radial grid integration of the TDSE 
using the matrix iteration method [8]. We employed a 
one-electron basis in an parametrized optimized effec- 
tive potential [9]. We used the linearly polarized XUV 
pulse £{t) = E g(t)smojt with the envelope g(t) rep- 
resented by the Nutall window function and centered 
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at t = 0. The following field parameters were chosen: 
E = 0.119 a.u. (corresponding to the peak intensity of 
5 x 10 14 W/cm 2 ), uj = 106 eV, T = 2tt/uj = 39 as and 
FWHM=182 as. Experimental field intensity was not re- 
ported by Schultze et al [5]. Given a typical high-order 
harmonic conversion efficiency of 10~ 6 [10], the presently 
chosen XUV field strength is most certainly larger than 
the one used experimentally. With this choice, our cal- 
culation is guaranteed to account for non-perturbative, 
with respect to the field, ionization effects if this effects 
were to be sizable. Other XUV field parameters used in 
the present work were identical to those used in the ex- 
periment. The XUV pulse described above is shown on 
the top panel of Fig. 1 (dotted black line). The pulse is 
truly off outside the interval ±Ti, where T\ ~ 5T which 
is about twice the FWHM. 

The solution of the TDSE satisfies the initial condi- 
tion <f(r,i < 7\) = <pi{r) which corresponds to a bound 
electron state on the atomic shell i to be ionized. So the 
shell index i is implicit in the following but omitted for 
brevity. The wavepacket representing the photoelectron 
ejected from a given shell is defined as 

^^EJ^WwiW^ dk , (1) 

where akh{t) — e lEkt {xkL\^{t)) are the projec- 
tion coefficients of the solution of the TDSE on 
the continuum spectrum of the atom. The con- 
tinuum state Xfei( r ) — Rki{r)Yi,(r / r) is the prod- 
uct of the radial orbital with the asymptotic Rki oc 
sin [kr + 5i(k) + l/fcln(2fcr) — Itt/2~\ and the spherical 
harmonic Y^r/r) with L = l,m. The projection coeffi- 
cients dki{t) cease to depend on time for \t\ > T\ when 
the driving XUV pulse is off. 

There are two convenient indicators of the evolution 
of the wave packet (1). One is the norm given by the 
integral Nit) = J2l I dk l«fci(0| 2 - Thi s 

norm is plotted 

on the top panel of Fig. 1 with the red solid and green 
dashed lines for the wave packets originated from the 2 s 
and 2p sub-shells, respectively. For better clarity, these 
curves are scaled and over-plotted on the electromagnetic 
pulse. 

The figure shows clearly that the evolution of the 2s 
and 2p wave packets starts and ends at the same time 
without any noticeable delay. This is further visualized in 
the inset where the variation of the norm [N (t) — N (Ti)] / 
N(Ti) is plotted on an expanded time scale near the 
driving pulse end. Indeed, the norm starts deviating 
from zero with the rise of the XUV pulse and reaches 
its asymptotic value once the interaction with the XUV 
pulse is over. 

Another marker of the wavepacket dynamics is the 
crest position, defined as a location of the global max- 
imum of the electron density. The latter quantity is 
truly informative only when the electron is outside the 
atom and the wavepacket is fully formed, having one well- 
defined global maximum. On the bottom panel of Fig. 
1, we show the crest position of the 2s and 2p wave pack- 
ets propagating in time. This figure can be viewed as a 
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FIG. 1: Top: the norm of the wave packets N(t) (scaled 
arbitrarily) emitted from the 2s and 2p sub-shells is plotted 
as a function of time with the red solid and green dashed lines, 
respectively. The XUV pulse is over-plotted with the black 
dotted line. In the inset, the norm variation [N(t) — N(Ti)]/ 
N(Ti) is shown on an expanded time scale near the pulse end. 
Bottom: the crest position of the 2s and 2p wave packets is 
shown with the same linestyles. The crest position after the 
pulse end is fitted with the straight line which corresponds to 
the free propagation. In the inset, extrapolation of the free 
propagation inside the atom is shown. 

more realistic version of a somewhat idealized and sim- 
plified graph presented in Fig. 1 of Schultze et al [5] . We 
see that evolution of the norm and the movement of the 
crest commence and cease at about the same time. 

The movement of the crest becomes almost linear when 
the norm reaches its asymptotic value and the wave 
packet is fully formed. Once fitted with the linear time 
dependence r = k(t — to) + tq for large times t > T\ 
(shown as a dotted straight line) and back propagated 
inside the atom, the 2s wavepacket seems to have an ear- 
lier start time to than that of the 2p wavepacket. This 
difference is magnified in the inset. It is about 6 as at the 
origin ro = and about 4 as at the distance ro ~ 1 au 
which corresponds to the size of the valence shell of the 
Ne atom. We see that at the origin ig s < and t p > 
are shifted to the opposite direction with respect to the 
peak of the driving XUV pulse which sets the start time 
of the photoionization process. Thus the seeming (or 
apparent according to Ref. [5]) "time zero" of the wave 
packet, which is inferred by the backward time propaga- 
tion, is different from the physical (or real) "time zero" 
t = 0. 

The origin of this shift is most clearly elucidated within 
the perturbation theory (PT) framework, which should 
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be applicable under the present field conditions of a single 
photon transition with 7 ^> 1 [11]. Under these condi- 
tions, 



/oo 
{ X kL\z\&)e l( - E "-^ T £(T) dr 
-OC 



(2) 



Here we extended the integration limits outside the pulse 
duration and wrote the dipole matrix element (xfeAl^l^j) 
in the length gauge. By separating the angular and radial 
integration, we can present this matrix element in the 
reduced form 



(XkL\z\(f>i) oc C l ^ l m .dx(k) 



(3) 



where C|™ ; . TO . is the Clebsch-Gordan coefficient, A = l,i 
and d\(k) is real. With this definition, we can write 
o-kL oc —i d\{k)£{Ek — e*), where the Fourier transform of 
the XUV field £(u) = e iUT 8{r) dr is real for a sym- 
metric pulse that we presently consider. We note that 
Eqs. (l)-(3) are equivalent to Eqs. (S5)-(S8) of Schultze 
et al [5] given in their supporting online material. 

To describe the motion of the wave packet (1), we 
apply the usual saddle-point method. For each I, the 
crest of the wave packet is moving at large times t>T\ 
quasi-classically along the trajectory which is given by 
the equation: 



(4) 



Since the logarithm is a slowly varying function which can 
be absorbed into a constant, Eq. (4) describes a straight 
line r = k(t-t ) + ro with t Q = dSi(k)/dE\ k= ^^- . Thus 
the relative time delay between various photoionization 
channels is determined primarily by the derivatives of the 
corresponding elastic scattering phases [12]. 

The scattering phases 5i(k) of the photoelectron mov- 
ing in the field of a singly charged Ne + ion are shown on 
the middle panel of Fig. 2. The photoelectron ejected 
from the 2s shell has only one value of the angular mo- 
mentum I = whereas the 2p photoelectron can acquire 
two angular momenta I = and 2. The phases in the 
s- and p- waves are shifted downwards by n and 7r/2, re- 
spectively, for better clarity. On the top panel of the 
same figure, we display the asymptotic projection coeffi- 
cients ai(fc) for I = 0, 1 and 2 and m = 0. The centers 
of the energy distribution of the /-projected coefficients 
(indicated by the vertical dotted lines on the top panel) 
define the position of the energy derivative of the corre- 
sponding scattering phases dS[(k)/dE (indicated 

fc=V2£fo 

by the straight lines on the middle panel). We see that 
the energy derivatives of the s- and p-phases are neg- 
ative whereas that of the d-phase is positive. This is 
so due to the presence of the occupied s and p states in 
the Ne + ion which disturbs otherwise monotonic increase 
with energy of the Coulomb phase (the Levinson-Seaton 
theorem [13]). Because the 2p — > kd transition is strongly 
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FIG. 2: Top: Expansion coefficients |a^;| 2 plotted versus the 
photoelectron energy Ek = fc 2 /2 which is expressed in eV. 



Middle: the HF scattering phases. Bottom: 
RPAE dipole matrix elements aig[D\(k)]. 
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dominant over the 2p — > ks one, as is seen from the corre- 
sponding projection coefficients on the top panel of Fig. 
2, it is the fi-phase that determines the shift of the ap- 
parent "time zero" of the 2p wave packet relative the 
physical "time zero" t — 0. This shift is positive fot the 
2p wave packet and negative for the 2s one, in accordance 
with our observation displayed in the inset of the bottom 
panel of Fig. 1. 

So far, we confined ourselves with an independent elec- 
tron approximation and calculated the dipole matrix el- 
ements d\(k) and the scattering phases 5i(k) in the HF 
approximation [14, 15]. It is well known, however, that 
many-electron correlation modifies strongly the dipole 
matrix elements. The full account for this effect can be 
taken within the RPAE model [6] by solving a set of cou- 
pled integral equations 



D x (k) =d x (k)+Y,J dpD u (p) xAp) U vX (j>, k) 



(5) 



Here Xvip) = (w — i? p — + ie) _1 is the Green's function 
and U u \(p,k) is the Coulomb interaction matrix. The 
one-electron HF basis corresponding to the field of the 
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singly charged Ne + ion accounts for direct photoelectron 
interaction with its parent shell. It is therefore the inter- 
shell Coulomb interaction with v ^ A that should only 
be included into Eq. (5). Since the Green's function is 
complex, the dipole matrix elements D\{k) acquire an 
additional phase which is plotted on the bottom panel of 
Fig. 2. 

The HF phase derivatives alone account for the appar- 
ent "time zero" shift between the 2s and 2p ionization 
AtQ S ~ 2p = 6.2 as. The RPAE correction adds an extra 
2.2 as. In total, this accounts for the apparent "time 
zero" shift At 2 Q s ~ 2p = 8.4 as. Both the HF and RPAE 
phases are smooth functions of the photoelectron energy 
and their averaging over the bandwidth of the XUV pulse 
does not change these numbers in a noticeable way. The 
analogous values reported by Schultze et al [5] for the in- 
dependent electron model and the correlation correction 
are 4.0 as and 2.4 as, respectively. Both sets of calcu- 
lations are quite close and well below the experimental 
value of 21 ± 5 as. 

One could argue that the complete account for many- 
electron correlation within the TDSE, rather than adding 
this correlation ad hoc, could modify the present result. 
This is, however, unlikely given the nature of the RPAE 
which is a direct generalization of the HF method in 
the presence of an oscillatory external electromagnetic 
field [16]. The only approximation taken when deriving 
Eq. (5) is that at any instant of time the atomic wave 
function is an anti-symmetric product of one-electron 
functions. It is quite a robust approximation under the 
field parameters considered above. 

We also evaluated the time delay of the wave packet 
relative to the XUV pulse in other noble gases. In He, the 
wave packet emitted from the Is shell is delayed by ~ 2 as 
relative to the center of the XUV pulse. This follows from 
the independent electron HF calculation which returns a 
positive derivative of the p phase shift as there is no occu- 
pied p orbital in the He + ion. It is also confirmed by the 
correlated convergent close-coupling (CCC) model which 
is known to produce benchmark photoionization results 
for He in the XUV range [17]. It is to be compared with 
5 as delay reported for He by Schultze et al [5] . In heav- 
ier noble gases, Ar and Kr, the difference of the HF p— 
and d— phase derivatives becomes smaller as occupation 
of the ionic orbitals increases in line with the Levinson- 
Seaton theorem. In Kr, the d— phase derivative becomes 



negative as the 3d orbital is occupied. Accordingly, the 
time delay between the wave packets emitted from the ns 
and np valence sub-shells is getting smaller. When the 
HF and RPAE phase derivatives are combined, it results 
in 5.8 as delay in Ar and nearly zero delay in Kr around 
the 100 eV photon energy mark. 

In conclusion, we examined various effects leading to 
the shift between the apparent "time zero" of the photo- 
electron wave packets emitted from the 2 s and 2p shells 
in neon relative to the center of the XUV pulse which 
sets the timing of the photoionization process. We found 
that this shift is primarily due to the energy derivative of 
the HF elastic scattering phase shifts which differs signif- 
icantly for various partial waves. The RPAE correction, 
which accounts for many-electron correlation, is rather 
small and cannot explain the profound difference between 
the theoretical and experimental time delay. 

The apparent "time zero" is only meaningful when the 
wave packet is detected at large distances from the atom 
as in attosecond streaking experiments. This apparent 
"time zero" has little to do with the real time when the 
atomic photoionization begins which is fully determined 
by the driving XUV pulse alone. In this sense, the at- 
tosecond streaking is not informative on the early stages 
of the photoionization process. However, this technique 
allows one to determine the energy derivative of the quan- 
tum phase of the dipole matrix element [18], thus facili- 
tating the so-called complete photoionization experiment 
[19] . This is particularly important in those targets where 
the many-electron correlation is significant. 

The full potential of the attosecond streaking technique 
and its successful application in atomic collision physics 
can only be realized if the current strong disagreement 
between theory and experiment in Ne is resolved. The 
present study was not able to do so. Our simulations and 
analytic arguments indicate that the XUV pulse alone 
cannot produce such a large time delay and it is the 
streaking IR field that is most likely responsible for this 
effect. 
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